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Abstract. We demonstrate that there is only one physical pro- 
cess required to explain the spectrum and the variability of the 
radio source at the dynamical center of our Galaxy, Sgr A*, in 
the frequency range from ssl to ps 1000 GHz, namely optically 
thin synchrotron radiation that is emitted from a population of 
relativistic electrons. We attribute the observed variability to 
variable energy input from an accretion disk around Sgr A* 
into the acceleration of the electrons. 
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1. Introduction 

Some years ago Lesch et al. (1988) showed that the nonther- 
mal radio spectra of the Galactic Centre, including Sgr A* and 
the extended components (often referred to as Bridge and Arc, 
Brown & Liszt 1984, Reich et al. 1988), can be explained as op- 
tically thin synchrotron radiation due to a monoenergetic elec- 
tron distribution. Then it was assumed that Sgr A* is the source 
of the energetic particles, which propagate into the extended 
components. Now, due to increasing resolution, the spectrum 
of Sgr A* itself can be traced up to 10 13 Hz. For frequencies up 
to at least 10 12 Hz Sgr A* exhibits the same inverted spectrum 
as the Bridge and the Arc, namely S v oc v 1 ^ ( v. frequency, S v : 
flux density) (see sect. 2 and references there). This is the typi- 
cal spectrum of electrons following a monoenergetic relativistic 
distribution and radiating optically thin synchrotron radiation. 
Thus, we can now make a much stronger point in favour of 
Sgr A* being a source of monoenergetic electrons, indeed. 

In the following first we discuss measurements of the spec- 
trum of Sgr A* and its time variability (Sect. 2). Then (Sect. 3) 



we give arguments which favour the existence of a monoener- 
getic electron distribution in the vicinity of an accreting black 
hole and calculate (in Sect. 4) the energetics of such a radiating 
system. In Sect. 5 we discuss the amplification of the required 
magnetic fields in an accretion disk around a central black hole. 
Finally (Sect. 6) we demonstrate its applicability to the Galactic 
Centre and to Sgr A* in particular. 

2. The spectrum of Sgr A* and its variability 

The spectrum of Sgr A* is described by Zylka et al. (1992) and 
by Mezger (1994). Its variabilty is well established. In figure 1 
we summarize the range of observed variabilty from 1 .5 GHz 
to 86 GHz (full lines), and individual measurements at higher 
frequencies up to 375 GHz (when given in the literature, bro- 
ken lines show the error bars of the individual measurement; 
otherwise an asterisk indicates the measured value). The obser- 
vations are taken from Backer (1982, 1994), Zylka & Mezger 
(1988), Zhao et al. (1989, 1991), Zhao & Goss (1993), Zylka et 
al. (1992), and Dent et al. (1993). 

Additionally, in figure 1 we give two best fit lines (linear 
in log S v vs. log v). The steeper one represents all observa- 
tions given in the figure. Its dependency is S v ~ v a 38 . For the 
shallower one we took into account only the frequencies for 
which a variation is well established (full lines). This ensures 
that we exclude the chance that a single observation at a certain 
frequency accidentally gives an extremal value within its range 
of possible variability. For the latter representation we get a 
dependency of S v 



,0.35 
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3. A monoenergetic electron distribution 

The particle momentum or energy distribution function is the 
result of the interplay of energy gain and energy loss processes. 
Acceleration due to diffusive shock waves, magnetohydrody- 
namical (MHD) waves or magnetic reconnection shift the parti- 
cles upwards in energy space, whereas radiation processes like 
synchrotron or inverse Compton scattering shift the particle to- 
wards lower energies. In a steady-state situation the resulting 



Fig. 1. The range of observed variabilty of Sgr A* from 1.5 GHz to 86 GHz (full lines), individual measurements at higher frequencies up 
to 375 GHz (when given in the literature, broken lines show the error bars of the individual measurement; otherwise an asterisk indicates the 
measured value). The observations are taken from Backer (1982, 1994), Zylka & Mezger (1988), Zhao et al. (1989, 1991), Zhao & Goss (1993), 
Zylka et al. (1992), and Dent et al. (1993). 



equilibrium distribution function of relativistic electrons in the 
momentum range p > p M /a is given by (Lesch and Schlick- 
eiser 1987) 

f(p)~p a exp(-p/p M )- (1) 

a denotes the ratio of the rate of shock wave acceleration (or 
due to magnetic reconnection) to the rate of acceleration by 
MHD waves, pm denotes that momentum where the timescale 
for diffusive shock wave acceleration equals the timescale for 
inverse Compton and synchrotron losses. Equation (1) holds 
provided the escape time of particles t is much longer than the 
acceleration time, i.e., for a (quasi-)stationary situation. In the 
case of magnetic reconnection the maximum particle energy is 
given by the product of the length of the reconnection sheet 
and the electric field parallel to the magnetic field lines (Lesch 
and Reich 1992; Lesch and Pohl 1992). Strictly speaking, a 
monoenergetic distribution at p = p M will be maintained if 



PGmniPM ) + PLoss (PM ) = . (2) 

The physical mechanism for the accumulation of particles 
around pm in the case of shock wave and MHD wave accelera- 
tion is evident: at momenta p < pm momentum gain dominates 
momentum loss, p G ain + Ploss > and the electrons are con- 
verted upwards in momentum space. At momenta p > pm the 
opposite is true pcain+PLoss < and the electrons are convected 
downwards in momentum space. Without momentum diffusion 
the equilibrium distribution would be singular at p = pm, be- 
cause there would be no physical process to remove particles 
from p = pm ■ Since we allowed momentum diffusion of parti- 
cles this (^-function behaviour is smoothed out and an isotropic 
monoenergetic distribution function is formed. 

This argumentation is also applicable to acceleration by 
magnetic reconnection. The approaching magnetic field lines 
confine the particles in a current sheet. These particles can only 
escape along the field lines, where the electric field accelerates 



much slower than the motion along the field lines. Thus, the 
escape time is longer than their acceleration time and results 
in a monoenergetic electron distribution function, i.e., as a rel- 
ativistic electron beam (Lesch and Reich 1992). Relativistic 
electron beams are unstable with respect to the excitation of 
plasma waves. The waves isotropize the beam electrons and 
the initially anisotropic distribution function develops into an 
isotropic monoenergetic one. This isotropisation happens on 
the timescale (Achatz et al. 1990) 



"pe 



1/2 



(3) 



5.6 • 10 4 ^n^ is the electron plasma frequency, n e is the 



number density of thermal electrons, n x is the number density 
of relativistic electrons and m p (m e ) is the proton (electron) 
mass. 

In the following we assume that the isotropisation time is 
shorter than the energy loss time due to radiation and we will 
check this assumption a posteriore. 

In both acceleration scenarios (shocks and MHD waves 
or magnetic reconnection) an isotropic monoenergetic electron 
distribution appears. Such a distribution produces an inverted 
optically thin synchrotron spectrum with a spectral index around 
+ 1/3 which is in good agreement with the observed time aver- 
aged spectrum of Sgr A* for frequencies up to 10 12 Hz. 



form the formula for synchrotron radiation. For a monoener- 
getic distribution function the typical emitted frequency corre- 
sponds to the maximum frequency ^ max at which the spectrum 
exponentially cuts off (see fig. 2). 



Fig. 2. The optically thin synchrotron spectrum of an ensemble of rel- 
ativistic monoenergetic electrons (full line). The broken line demon- 



strates the quality of the approximation with S v 
below the cut-off. 



for frequencies 



4. The energetics of a radiating plasma with monoenergetic 
relativistic electrons 

In the following, we discuss the energetics of a radiating plasma 
with monoenergetic relativistic electrons, and use the observed 
luminosity L of Sgr A* of a few times 10 4 L Q (Mezger 1994, 
and references therein). The luminosity L of an object with 
emission volume V filled with radiating monoenergetic rela- 
tivistic electrons is given by (Lang 1978) 



L 



P( 7 )n t V, (4) 
where 

P(j) =1.6-10- 15 B 2 7 2 erg s" 1 (4a) 

denotes the emitted power of one electron with the Lorentz 
factor 7 emitting synchrotron radiation in a magnetic field B. 

We estimate the number density of relativistic electrons by 
assuming equipartition between energy density of relativistic 
electrons and magnetic energy density: 

B 2 
57T7m e c z 

The emission volume is taken from the 86 GHz observations of 
Krichbaum et al. (1993) who observed a source size of about 
0.34 milliarcseconds, which correspond to a linear scale of 
4.42 • 10 13 cm at a distance of the galactic centre of 8.5 kpc. 
If the emission region is an elongated structure, as indicated 
by the VLBI measurements (for details of the structure, see 
Krichbaum et al. 1993) the volume is V ~ 2.7 • 10 41 cm 3 . 



^max — 7 ^ce; 



(6) 



where the electron gyrofrequency is given by v ce = 5^77 • Thus, 
the particle Lorentz factor is 
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Inserting (4a), (5), and (7) into 4 we obtain 

L^4.65-10- 14 ^5 7/2 ^ 
which results in a magnetic field strength 

B ~ 13. 8G 



(7) 



(8) 
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With Eq. (9) we obtain from Eq. (7) the typical Lorentz factor 
of the radiating electrons 
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which corresponds to an energy of about 250 MeV. 
The energy loss time of these particles is given by 



(10) 
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Next, we check whether the isotropisation time [Eq. (3)] is 
shorter than ^loss as we assumed in the beginning. According to 
Eq. (3) t lso ~ 7.5 • 10~ 4 sec n l J 2 /n r . Since Sgr A* is embedded 
in a ffll-region, named Sgr A West (Lacy et al. 1980; 1982) 
we choose as an example the density of a compact HH-region 
n e ~ 10 6 cm -3 . We obtain a fraction of a second fort; so , which 
is much shorter than any reasonable energy loss time [eq. (11)]. 
Hence our assumption of an isotropic monoenergetic electron 
distribution is exceedingly likely. 

5. Accretion disk and magnetic field amplification 

Finally, we have to check whether our result for the magnetic 
field strength is reasonable in terms of an accretion disk model, 
including a massive black hole of the order of 1O 6 M in its 
center. 

Magnetic fields have long been considered an important 
element in the dynamics of accretion disks, primarly as a mech- 
anism for supplying internal stresses required for efficient angu- 
lar momentum transfer (Eardley and Lightman 1975; Ichimaru 
1977). According to these models magnetic fields are amplified 
within the inner portion of an accretion disk by the joint action 
of thermal convection and differential rotation along Keplerian 
orbits. Field amplification will then be limited by nonlinear 
effects; as a consequence of buoyancy, magnetic flux will be 
expelled from the disk, leading to an accretion disk corona con- 
sisting of many magnetic loops where the energy is stored and 
probably transferred via magnetic reconnection into heat and 
particle acceleration, leading to a state in which loops contain 
very hot, relatively low-density plasma (compared to the disk). 

Magnetic field amplification due to differential motion of a 
conductive disk medium is usually described in the MHD limit 
through the induction equation for the magnetic field (cf . Parker 
1979) 



SB „ 



[V 



B] + V 2 B, 

47T<7 



(13) 



where B is the magnetic field induction, v the velocity of the 
fluid, and a the electrical conductivity. Due to the high tem- 
peratures in the disk (up to 4 • 10 4 K; Zylka et al. 1992) the 
electrical conductivity is high. Therefore, magnetic field evolu- 
tion via the plasma differential motion is well described by Eq. 
(13) even if we neglect the classical (collisional) dissipation of 
the magnetic field implied by the last term on the right-hand 
side of Eq. (13). In the case of Keplerian differential rotation, 
Eq. (13) describes (assuming an axisymmetric geometry) the 
amplification of a toroidal magnetic field Bj in the presence of 
a seed poloidal magnetic field B P : 



dB T _ dCl 
Hf- R dR Bp ' 



(14) 



material about the black hole with mass M BH . 

The buildup of magnetic fields within the disk is limited 
by nonlinear effects related to convection. Since convection 
takes place primarly perpendicular to the plane of the disk, 
we shall assume differential rotation to remain the dominant 
mechanism for the toroidal magnetic field Bj amplification; the 
amplification of the poloidal field Bp will then be dominated by 
convection-mediated effects. For convection cells whose aspect 
ratio is ~ 1, the poloidal magnetic field spatial scale will then 
be of the order of the convective cell size H , equal to the half- 
thickness of the disk. To describe the amplification of Bp we 
invoke magnetic flux conservation: 



Bp — — Bt. 

R 



(15) 



The maximal toroidal field strength is given by the energy 
density of the turbulence in the disk (Galeev et al. 1979) 



r2 

n T 2 
— ~ n e m p c s , 
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(16) 



where c s = . / ^f- denotes the sound velocity. 

V P 

The effective temperature in the accretion disk at a radius 
R from the centre of the accreting black hole T accr (i?) is given 
by (see, e.g., Frank, King & Raine 1992) 



Tracer {R) — Tin 
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(17) 



with the temperature at the disk's inner edge R in ~ 3R$ of 



Tin = 7.1 • 10 4 K 
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(18) 



A temperature of 71,000 K corresponds to a sound velocity of 
2.4- 10 6 cms" 1 . 

For eqs. (17) and (18), we use the following approxima- 
tions: the disk is stationary (i.e., the radial mass flow rate M is 
constant); the non-relativistic limit of disk physics applies; and, 
the disk's inner edge can be described with the same physics as 
regions within the disk (i.e., we neglect the term coming from 
the radially inner boundary condition). It is known that the latter 
two approximations introduce only errors that are limited to the 
radially innermost disk regions and that correspond to factors 
of the order of 1, i.e., these approximations will not influence 
our conclusions. 

With the same parameters we obtain for the density a central 
value of about n e ~ 4 • 10 cm" - 3 (R/R in )- 3 / 2 . Altogether we 
get from Eq. (18) and (16) the maximum toroidal field strength 
of about B T ~ 30 G or a poloidal field of about Bp ~ 1 G, 
respectively. 



Since the only available radio observation which can be used 
for innermost parts of Sgr A* is the VLBI measurement by 
Krichbaum et al. (1993) we have to deal with an elongated 
structure, which might be interpreted as ajet (Falcke et al. 1993). 
Hence, we will take the poloidal field as the field strength which 
make the relativistic electrons radiate synchrotron emission. A 
value of a few Gauss is very close to the value we obtained from 
Eq. (9). 

Averaged over time, the spectrum of Sgr A* is very close 
to a S v ~ v x l 3 law (Fig. 1 in sect. 2). But it is important to 
note that individual spectra vary on timescale of weeks (Zhao 
et al. 1989, 1991; Zhao & Goss 1993) and sometimes differ 
considerably from a S v oc v 1 ^ dependency. 

Together with the above shown consistency of the observed 
and the theoretical magnetic field, this leads us to the conclusion 
that what we observe in Sgr A* is indeed driven by a monoen- 
ergetic distribution of electrons emitting synchrotron radiation. 
The variability is then due to a variable energy input from the 
accretion disk. As shown above [eq. (11)] the energy loss time 
is of the order of a few days. Within such a time scale, the 
electron population can react upon a variable supply of energy 
for the acceleration process. 

Thus, the fact that the time averaged spectrum indeed shows 
a!/ 1 ' 3 dependency tells us that the timescale within which the 
energy input into the acceleration process may vary must be of 
the same order as the loss time scale. 

Were the timescale of the input very much shorter then the 
electron distribution would "see" only the averaged energy input 
rate, i.e., we would not expect any variations of the spectrum 
at all. This argument again relys on an accretion rate that is 
constant on the average. It also takes not into account other 
reasons for variability (e.g., refractive interstellar scintillation; 
this process could account for [some] variability at the lower 
frequencies [< 10 GHz], but hardly for the higher ones; see, 
e.g., Zhao et al. 1990). 

Were the time scale very much longer then again we would 
not expect deviations from the proportionalty to v x l 3 as then the 
electrons could follow the changes in a quasi-stationary fashion. 

With the above estimate for the sound velocity, we find that 
characteristic viscosity time scales in the central regions of an 
a accretion disk (i.e., at radial distances of a few to several 
Rs) are of the order of a few days, i.e., of the same order as 
the energy loss time scale of the electrons. This demonstrates 
that the observed time scale of spectral variations and both 
the theoretically deduced time scales (electron energy loss and 
viscous changes in a disk) are in good agreement with each 
other. 



7. Conclusions 

We have shown that the spectrum of Sgr A* in the frequency 
range of log(^/Hz) ps 9 ... 12 and its variability can be inter- 
preted as being due to a single underlying physical mechanism, 



that emit synchrotron radiation. 

For our arguments the variability of the spectrum is crucial. 
Only the near coincidence of the two relevant time scales (ac- 
cretion time scale and energy loss time scale of the electrons) 
allows a variability in which individual spectra may very well 
deviate from the v 1 ^ law while - on the average - the source 
knows about the steady state solution. 

In a next step one has to look into the details of the time 
dependent evolution of such an electron distribution. This then 
will also allow to decide whether only the high frequency vari- 
ations (> 10 Ghz) are due to intrinsic source variations. For 
lower frequencies, in our model, variability could also be due 
to source variability. On the other hand, Zhao & Goss (1993) 
give arguments that favour scattering of the radiation by the 
ISM along the line of sight as the reason for variability at the 
lower frequencies. It seems as if this can be clarified only by a 
time dependent model. But such a model is beyond the scope 
of our present work that is presented in this paper. 

Supporting evidence for Sgr A* being a source of monoen- 
ergetic relativistic electrons comes from the fact that also the 
Bridge and the Arc in the Galactic Center show the same kind 
of spectrum (Lesch et al. 1988). 
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